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SUMMARY 


One  problem  with  using  the  MICRAD  Attitude  Reference  System 
(MARS)  to  provide  error  signals  to  the  Vertical  Seeking  Seat  (VSS) 
control  system  is  that  the  MARS  outputs  vary  significantly  depending 
on  whether  the  local  earth's  surface  is  land  or  water.  This 
variation  changes  the  control  system  gain  resulting  in  degraded 
performance.  As  a  result,  efforts  were  initiated  to  develop  alter¬ 
natives  and/or  improvements  to  control  laws  studied  previously  at 
the  Naval  Weapons  Center  (NWC).  An  alternate  control  law  was 
developed  and  evaluated  using  a  digital  simulation.  Results 
indicated  that  this  control  law  provided  good  performance  and 
adequately  compensated  for  variations  in  the  MARS  outputs. 
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INTRODUCTION 


The  Maximum  Performance  Ejection  Seat  (MPES)  or  Vertical  Seeking 
Seat  (VSS)  escape  system  requires  a  vertical  reference  to  provide 
commands  for  the  control  system.  Previous  tests  of  the  MPES/VSS 
system  were  performed  using  a  strap-down  vertical  reference  system 
which  was  initialized  at  the  planned  "launch"  attitude.  Although 
this  attitude  hand-off  is  satisfactory  for  initial  tests,  the 
operational  system  should  use  an  autonomous  reference  system.  One 
such  system  which  is  to  be  used  in  future  flight  tests  ts  the  Hlcro- 
wave  Radiometric  (MICRAD)  Attitude  Reference  System  (MARS). 

The  NARS  contains  a  reference  load  and  four  antennas  which 
measure  the  natural  radiation  emanating  from  objects  within  their 
field  of  view,  such  as  the  earth  or  sky  or  combinations  of  the  two. 
The  outputs  of  the  MARS  are  an  indication  of  the  relative  tempera¬ 
tures  of  the  load  and  those  sensed  by  each  antenna.  These  four 
measures  of  relative  temperature  are  to  be  processed  by  control  laws 
whose  outputs  are  used  to  drive  angular  velocity  feedback  control 
loops.  One  area  of  concern  Is  that  the  significant  difference 
between  the  microwave  temperatures  of  land  and  water  may  change  the 
gain  of  the  control  laws,  thus  making  the  system  response  dependent 
upon  the  local  make-up  of  the  earth's  surface. 

The  objective  of  this  effort  is  to  develop  alternatives  and 
improvements  to  control  laws  studied  previously  at  the  Naval  Weapons 
Center  (NWC) .  Specifically,  the  problem  of  compensating  for  the 
different  temperatures  of  land  and  water  wi>'  te*  a'Jdressed.  The 
resulting  control  laws  are  to  be  suitable  for  implementation  in  r»>e 
on-board  microprocessor.  Although  the  candidate  control  laws  need 
not  be  optimized  for  all  flight  conditions,  they  will  be  evaluated 
for  several  flight  conditions  using  a  six  Degree  of  Freedom  program 
provided  by  NWC. 
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CANDIDATE  CONTROL  LAWS 

Three  control. laws  have  previously  been  studied  at  NWC  for  use 
with  the  MARS.  These  laws  and  their  origin  are  briefly  discussed  here 
for  completeness.  In  the  case  of  the  Heuristic  Control  Law,  an 
attempted  derivation  yielded  a  modified  control  law  which  is  also 
considered  as  a  candidate. 
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DIRECTION  COSINE  CONTROL  LAW 

The  Direction  Cosine  Control  Lew  (DCCL)  is  based  on  commanding 
body  angular  velocities,  q  and  p,  proportional  to  the  cosines  of  the 
angles  between  the  nadir  and  the  x  and  y  body  axes  respectively.  One 
method  of  viewing  and  implementing  the  DCCL  with  the  MARS  outputs  is 
shown  here.  Consider  the  component  of  a  vector,  representing  the 
temperature  difference  between  the  nadir  and  zenith,  in  the  seat  body 
x  -  y  plane  as  shown  in  Figure  1.  When  the  seat  is  aligned  with  the 
vertical,  the  vector  component  AT  is  zero.  This  suggests  that  the 
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FIGURE  I.  Seat  Body  x  -  y  Plane  Showing  Temperature  Difference  and 
Desired  Angular  Velocity. 


scat  be  rotated  about  an  axis  perpendicular  to  AT  in  the  x  -  y  plane. 
Thus  the  commanded  angular  velocity,  shown  in  Figure  I,  can  be 
expressed  in  terms  of  AT. 


If 


then 


fiT  o  q  x  •  p  y 
W  “  p  x  +  q  y 


(1) 
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In  order  to  obtain  p  and  q  as  functions  of  the  antenna  outputs,  we 
must  determine  components  of  AT  along  the  antenna  axes. 


a£  =  -  q  2 

a^  =  (q  cos  60°  "  p  cos  30°)  3  (2) 

a^  e»  (q  cos  60*  +  p  cos  309)  4 


These  equations  can  be  rearranged  to  yield 
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a4  1  a3 

/r 
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(3) 


which  define  the  corananded  angular  velocity  in  body  axes  in  terms  of 
the  normalized  antenna  outputs  3j.  It  was  shown  in  Reference  1  that 
if  the  normalized  antenna  outputs  a2,  a^,  and  a^  are  cosines  of  the 
angles  between  the  respective  antenna  and  the  nadir,  then  Equation 
(3)  is  equivalent  to  the  OCCL  if  p  and  q  are  the  commanded  roll  and 
pitch  rates  respectively.  In  general  this  will  not_  be  the  case 
since  the  antenna  outputs  will  not  be  normalized  to  a  cosine  function.. 
Therefore,  Equation  (3)  will  be  referred  to  as  the  KARS  Direction  Cosine 
Control  Law  (HLCu*.) . 


HEURISTIC  AND  MOOIFIEO  HEURISTIC  CONTROL  LAW 
The  Heuristic  Control  Law  (HCL)  is 
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a  a4  '  °3 
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where  p  and  q  are  the  commanded  angular  velocities  and  the  antenna 
outputs  a2,  a3,  and  a^  are  considered  to  be  normalized.  An  expression 
very  similar  to  Equation  (4)  can  be  obtained  by  resolving  the  antenna 
outputs  into  the  x  and  y  body  axes. 


From  Figure  1 
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Now  define  che  cotwaanded  angular  velocity  as 
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Note  that  Equations  (U)  and  (5)  differ  by  a  constant  factor  in  the 
expression  for  p.  Therefore,  Equation  (5)  will  also  be  considered 
as  a  candidate  control  law  and  will  be  referred  to  as  the  Hodified 
Heuristic  Control  Law  (HHCL). 


PERHUTATION  CONTROL  LAM 

The  Permutation,  or  table  look-up.  Control  Law  (PCL)  is  based 
on  sorting  the  four  KARS  antenna  outputs  in  order  of  coldest  to 
warmest  and  using  predetermined  angular  velocity  commands  for  each 
perrautetion  of  the  antenna  outputs.  Since  the  resolution  of  the 
resulting  on-off  control  logic  is  significantly  less  than  the 
desired  system  accuracy,  additional  proportional  control  logic  is 
required  when  the  seat  is  near  its  desired  attitude.  A  detailed 
description  of  the  PCL  can  be  found  in  Reference  1. 


analysis  of  control  laws 


Each  of  the  previously  described  control  laws  is  briefly 
analyzed.  Initially,  the  proportional  control  laws,  that  is,  the 
Direction  Cosine,  Heuristic,  and  Hodified  Heuristic  Laws  are 
evaluated  with  regard  to  their  characteristics  using  un-normal ized 
antenna  outputs  with  varying  nadir  temperatures.  Next,  the  error 
between  the  consnanded  angular  velocity  axis  defined  by  each  of  the 
control  laws  is  compared  with  an  axis  which  minimizes  the  rotation 
angle  required  to  align  the  seat  with  the  vertical.  Finally,  a 
similar  comparison  of  the  permutation  law  commands  with  those 
required  for  a  minimal  rotation  axis  is  made.  Based  on  this  analysis 
a  control  law  is  selected  for  further  evaluation  using  the  NWC  six 
Degree  of  Freedom  (DOF)  digital  simulation. 

EFFECTS  OF  MARS  OUTPUT  CHARACTERISTICS 

As  discussed  previously,  the  MARS  antennas  measure  temperature 

and,  therefore,  their  outputs  are  significantly  dependent;  on  whether 
the  sensor  is  over  land  or  water.  This  causes  large  gain  or  scale 
factor  variations  in  the  proportional  control  laws,  if  un-normalized 
antenna  outputs  are  used  directly,  which  may  degrade  the  system 
response.  One  solution  to  this  problem  is  to  normalize  the  antenna 
outputs.  However,  this  may  be  difficult  since  the  nadir  and  zenith 
temperatures  are  not  kno*»m  a  priori  and  may  not  both  be  measured, 
depending  on  the  seat  trajectory.  Insight  is  obtained  by  evaluating 
the  control  laws  using  an  expression  for  the  normalized  antenna 
outputs  in  terms  of  the  actual  outputs.  Assume  that  the  actual 
output  of  an  antenna  at  the  horizon  A  is  the  average  of  the  actual 
outputs  of  antennas  pointing  at  the  nadir  Ar  and  zenith  A^. 


How,  assuming  WARS  outputs  high  (low)  voltages  for  low  (high)  temp¬ 
eratures.  define  the  normalized  output  a-  of  antenna  i  as 


Where  Aj  is  the  actua'i  output  of  antenna  i.  Substituting  Equation 
(6)  into  (7)  and  rearranging  yields 


2A,  -  <A„  ♦  V 


(8) 


The  normalized  antenna  outputs  defined  by  Equation  (8)  are  used 
directly  in  the  HDCCL  and  HCL,  that  is,  in  Equations  (3)  and  (4) 
respectively.  The  resulting  HDCCL  is 

2(AJj  -  A}) 


<\  -  V  * 


r(A,  ♦  A,  -  A.)  -  (An  *  a2) 


2* 

2<Az  -  \> 


(9) 


and  the  HCL  is 


2  tA4  "  V 

p  “ 

z  n 
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q  *  -r-nr  ( V 

z  n 


(10) 


The  HHCl  is  not  listed  as  it  is  very  similar  to  the  HCL.  Mote  that 
in  Equations  (9)  and  (10).  A2,  *y  and  are  the  un-normalized 
antenna  outputs  and  therefore  these  equations  wilt  be  referred  to 
as  un-normal Szed  control  laws, 

A  comparison  of  Equations  (4)  and  (10)  indicates  the  normalized 
and  un-normal i zed  HCL’s  are  related  by  a  factor  of  2/(A2  -  A^) .  The 
relationship  of  the  corresponding  KDCCl’s,  Equations  (})  and  (9),  is 
sorHSwhat  more  complicated,  implementation  of  Equations  (9)  and  (10) 
requires  knowledge  of  (A,  -  Zj  and  the  former  also  requires  know¬ 
ledge  o'  (A^  ♦  A^).  As  mentioned  previously,  the  exact  values  of 

A  and  A  are  not  known  a  priori.  However,  it  may  be  possible  to 
n  2 

estimate  them  using  the  fact  that  is  relatively  constant  under 
varying  am»tcoro logical  conditions.  Thus  the  use  of  a  predetermined 


constant  to  approximate  Az  may  be  feasible.  Aiso  helpful  is  the 
following  approximation  for  the  actual  output  of  an  antenna  pointing 
at  the  horizon. 


A  = 


A2  +  A3  +  A^ 
3 


(11) 


Equation  (6)  can  be  rearranged  to  yield 

A  ♦  A  3  2 
n  z 


(12) 


where  A  is  obtained  using  Equation  (11).  Now  Equation  (12)  can  be 
used  along  with  a  predetermined  estimate  of  designated  as  A^,  to 
provide 

Az  -  An  -  2  (Az  -  A)  (13) 


where  A  is  again  obtained  using  Equation  (11).  Equations  (11). 
(12),  and  (13)  can  be  used  to  estimate  the  quantities  required  for 
implementation  of  the  un-'noraal ized  control  laws  as  follows.  The 
HDCCL  1 S  now 


and  the  HCl  is 


<*(,  -  y 
-  »)  fT 

(Aj  ♦  A^  *  ftj) 

2 (A,  -  A) 


04) 


A4  '  *3 

w  J  -  ■ 

A  -  A 
z 


05) 


A  -  A 
z 


VA 


-  V 


where  A  is  defined  by  Equation  (51).  The  KHCl  is  the  saa*e  as 
Equation  (IS)  except  the  fonser  has  an  additional  factor  of  .'3/2  in 
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the  expression  for  p.  The  effectiveness  of  the  above  approximations, 
Equations  (it),  (12),  and  (13),  in  dealing  with  varying  nadir 
temperatures  and  consequently  varying  scale  factors  will  be  eval¬ 
uated  using  scaled  MARS  data  in  the  NWC  6-DOF  simulation. 

The  PCI  will  require  the  same  type  of  compensation  derived 
above  since  it  requires  the  use  of  a  proportional  control  law  for 
small  offsets  from  the  vertical.  For  large  offsets  no  such  compen¬ 
sation  is  required  as  the  PCL  uses  the  relative  order  of  the 
antenna  output  magnitudes  to  determine  the  angular  velocity  commands. 


DIRECTION  OF  COMMANDED  ANGULAR  VELOCITY  VECTOR 

One  interesting  characteristic  of  the  control  Jaws  is  the 
difference  between  the  directions  of  the  commanded  angular  velocity 
vector  and  the  vector  which  minimizes  the  angular  rotation  required 
to  align  the  seat  with  the  vertical,  assuming  the  body  angular 
velocity  r  e  0.  The  direction  of  the  Minimum  Rotation  Axis  (MRA) 
is  found  as  the  cross  product  of  body  z  axis  and  the  inertial  z 
axis  expressed  in  body  axes.  Using  the  standard  seat  coordinate 
systems  and  Euler  angles. 


0 

ro' 

MRA  ■ 

0 

1 

— 

X  CB/I 

0 

J. 

where  is  the  Direction  Cosine  Matrix  (DCM)  between  body  and 

inertial  axes.  This  yields 


sin  $  cos  0’ 

"C23 

MRA  - 

-  sin  0 

* 

0 

*■  •* 

.0 

(16) 


where  0  and  $  are  the  pitch  and  roll  Euler  angles  of  the  seat,  and 

C?2  ond  C22  «re  components  of  the  DCM.  Note  that,  by  definition, 

C,_  and  C^_  are  the  direction  cosines  between  the  inertial  z  axis 

13  o 

and  the  body  x  and  y  axes  and  are  therefore  equal  to  the  outputs 

q  and  p  respectively  of  the  ideal  DCCL.  In  order  to  evaluate  the 

control  law  outputs,  expressions  from  Reference  2  which  provide  the 
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angle  between  each  antenna  axis  and  the  zenith  were  used.  These 
expressions  are  repeated  here  for  convenience. 

Zj  =>  cos  1  (cos  $  cos  0) 

z2  =  90°  +  0  (17) 

=  cos  1  [(sin  0-3  sin  $  cos  e)/2] 

=  cos  1  [(sin  6  +  3  sin  41  cos  e)/2] 

where  z;  is  the  angle  between  antenna  i  and  the  zenith.  The  actual 
antenna  outputs  are  obtained  by  using  the  above  angles  and  inter¬ 
polating  the  MARS  output  data  for  Antenna  #1  in  Appendix  A.  A  list 
of  the  antenna  outputs  for  several  seat  orientations  is  given  in 
Table  1.  A  comparison  of  the  corresponding  control  law  outputs 
and  MRA  is  presented  in  Table  2,  where  the  orientations  of  the 
commande  .  ar  velocity  vector  and  the  MRA  are  the  angles  from 
the  x  body  axis  of  the  seat.  That  is 

/.MRA  •  tan"1  (|?~-&)  (18) 

from  Equation  (16),  and 

L.  Wc  -  tan"1  (£)  09) 

from  Figure  1. 

It  should  be  noted  that  the  data  in  Table  2  were  obtained  using 

the  idealized  control  laws  in  Equations  (S)  and  (10).  None  of  the 

previously  discussed  approximations  (Equations  (11),  (12),  and  (13)) 

were  used,  as  the  acutal  values  of  A  and  A„  were  obtained  from 

n  2 

Appendix  A.  Since  these  quantities  cancel  in  Equation  (19)  for  the 
HCL  and  the  HHCL,  the  use  of  approximations  would  only  affect  the 
data  from  the  HOCCL.  The  PCI  was  supplemented  with  the  ideal ized 
DCCL  to  provide  proportional  control  in  the  appropriate  regions. 


This  DCCL,  which  yields  the  MRA,  was  used  so  that  the  PCL  could 
be  evaluated  independently  of  errors  produced  by  the  other  proport¬ 
ional  control  Jaws.  Evaluation  of  the  data  In  Table  2  indicates 
that  the  control  laws  can  be  ranked  in  terms  of  the  difference 
between  their  rotation  axis  and  the  MRA  as  follows;  closest  is 
the  MHCL,  followed  by  the  MDCCL,  HCL,  and  PCL.  Although  the  small 
differences  between  the  first  three  rotation  axes  are  not 
presently  considered  critical,  their  importance  may  increase  as 
the  system  performance  is  optimized  u r  if  three  axis  control  is 
implemented.  Also  recall  that  the  MDCCL  axis  will  be  altered  by 
the  approximations  required  to  implement  the  un-normal tzed  control 
laws.  For  these  reasons  the  following  MHCL  was  selected  for 
evaluation  using  the  NVC  6-OOF  simulation. 


where  the  subscript  "c"  has  been  added  to  indicate  commanded  values 
and  K  is  a  constant  gain. 
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ADDITIONAL  COMPENSATION 


Previously  the  proportional  control  laws  were  evaluated  by 
comparing  the  directions  of  the  commanded  angular  velocity  with  a 
desired  direction  assuming  the  body  angular  velocity  r  was  held 
to  zero.  In  actuality,  these  control  laws  may  cause  limiting  in 
some  component  of  the  system,  thus  losing  their  proportional  char¬ 
acteristics.  Also,  r  will  not  be  2ero,  Methods  of  compensating 
for  these  effects  are  discussed  below. 

ANGULAR  VELOCITY  LIMITS 

The  seat  system  contains  two  independent  angular  velocity  feed 
back  control  loops  which  are  driven  by  the  control  law  outputs  p 

v 

and  qc.  If  any  portion  of  this  system  reaches  an  operating  limit, 
the  proportional  characteristics  will  be  lost,  resulting  in 
directional  divergence  of  the  angular  velocity  vectors  provided  by 
the  control  law  and  achieved  by  the  system.  One  method  of  prevent¬ 
ing  this  is  to  limit  the  maximum  values  of  the  angular  velocity 
commands  p  and  q  .  Care  should  be  taken  in  implementing  these 
limits  so  as  to  preserve  the  direction  of  the  commanded  angular 
veioci ty. 

Consider  the  situation  where  the  angular  velocity  commands 

provided  by  the  control  law  exceed  the  desired  limits  as  shown  in 

Figure  2.  In  this  case  the  control  law  output  W  exceeds  the  limit 

for  q  .  If  q  is  simply  limited  to  W.  and  p  Is  unchanged,  then 

the  axis  of  the  commanded  angular  velocity  will  have  changed  as 

indicated  by  AW  in  the  figure.  Another  approach  Is  to  preserve  the 

angular  velocity  axis  by  changing  the  command  from  W  to  U  1 .  This 

can  be  accomplished  as  follows,  assuming  Is  the  desired  limit  on 

the  control  law  outputs  p  and  q  . 

c  c 

If  q  >  W.  >  p  or  q  >  p  >  U 

C  m  C  c  c  * 
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But  if 


then  let 


Pc  >  WA  >  \ 


Pc'  *  W* 


qc 

)  '  »  — 
c  p„ 


qc  <  Wjt  >  Pc 


Pc  “  pc 


°r  Pc  >  qc  >  W4 


V  •  qc 

where  p  1  and  q  '  are  the  new  angular  velocity  commands  after  limiting. 

c  c 


Ut  «  Desired  limits 


W  «  Control  law  Output 


P  Wfl 

*c  t 


i  -  *  «  *  • 


\-V 

M.V' 


O  >  V 


r  •  **  ’/ 
V  V.V 

i 

l."  *  m 

■■ 


W* 

►» 

fV-V' 

l\v‘> 

*’>;.■ 

SVAn 

r 


FIGURE  2,  Angular  Velocity  limits 


PHASE  LAG 


The  present  system  provides  for  control  of  two  of  the  seat 

body  axis  angular  velocities  p  and  q.  Although  the  third  angular 

velocity  r  Is  not  controlled,  It  is  measured  and  may  be  used  to 

provide  additional  compensation.  The  effects  of  r  on  the  commanded 

values  of  p  and  q  can  be  visualized  using  figure  3 •  The  outputs  of 

the  control  law  p  and  q  are  used  to  form  the  commanded  angular 
c  c 

velocity  vector  W  In  the  x  -  y  plane  of  the  seat  body.  If  r  is 
c 

zero  and  the  control  system  provides  an  angular  velocity  W  pro¬ 
portional  to  W.,  then  the  direction  of  W  will  be  constant  and  its 
c  c 

magnitude  will  diminish  as  the  seat  approaches  its  desired 

attitude.  In  this  case,  W  and  W  are  essentially  col  Inear  if  the 

c 

gain  and  phase  characteristics  of  the  two  angular  velocity  feed¬ 
back  control  loops  are  similar,  at  the  frequencies  of  interest, 
causing  p  and  q  to  be  in  phase.  This  is  a  reasonable  assumption 
as  these  control  loops  are  similar. 


FIGURE  3>  Relationship  of  Angular  Velocity  and  Its  Commanded 
Value  with  r  *  0. 

However,  If  r  is  not  zero,  Wc  will  rotate  j_n  the  x  -  y  plane  at  an 
angular  frequency  which  is  equal  In  magnitude  but  in  the  opposite 
direction  of  r.  This  rotation  Introduces  a  phase  shift  between  p 

V 

and  q  which  causes  the  U  and  U  axes  to  separate  as  shown  in 
c  c 

Figure  *«.  Now  U  lags  W£  by  the  phase  angle  e  as  both  vectors 
rotate  in  the  plane  of  the  figure.  If  the  seat  angular  velocities 
p,  q  and  r  are  well  below  the  bandwidth  of  angular  velocity  control 
loops,  then  the  phase  angle  e  may  be  reasonably  approximated  by  the 


nwdlfied  W  by  e  It,  W,  will  be  aligned  with  Its  desired  direction. 
This  is  accomplished  using  the  following  transformation 


-e  1  0 

0  0  1 


where  *  indicates  the  modified  commands.  Substituting  Equation  (25) 
into  (26)  yields 


pc*  ■  Pc  +  *05  r  qc 


qc*  ■  qc  '  *05  r  Pc 
where  r  is  in  radians/ sec. 

It  should  be  noted  that,  since  this  compensation  was  derived 
for  steady  state  conditions,  it  will  not  be  exact  under  other 
operating  conditions. 


EVALUATION  OF  CONTROL  LAWS 


The  HHCL  was  selected  for  further  evaluation  using  the  NWC 
6- DOF  digital  simulation..  Since  the  scope  of  this  effort  did  not 
Tiliow  for  optimization  of  the  HHCL,  some  other  method  was  required 
for  determining  the  gain  K  in  Equation  (20).  The  present  angular 
velocity  feedback  control  loops  were  designed  assuming  they  were 
driven  by  the  ideal  OCCL.  Therefore,  a  logical  approach  is  to 
match  the  gain  of  the  HHCL  and  the  Ideal  OCCL,  for  small  pertur¬ 
bations  from  vertical,  by  proper  selection  of  K.  This  Is  accom¬ 
plished  with  the  aid  of  the  block  diagram  in  Figure  5*  For  the 
ideal  DCCL 

K  «  1 

and  the  effective  gain  Is 

P  _  0.5**  _  6  rad/sec  / 

K“ ’ 'm  “Fid—  (28) 

The  corresponding  gain  for  the  HHCL,  calculated  using  Equation  (20) 
and  the  MARS  data  in  Appendix  C  for  wannest  nadir,  is 

K(2.70)|^-  K(2.70>^  U9) 

Equating  Equations  (28)  and  (29)  yields 

K  -  0.37  (30) 

This  uaiue  of  K  was  used  in  Equation  (20)  and  implemented  in  the 
simulation  as  shown  in  Figure  5.  The  response  of  the  resulting 
HHCL  was  compared  with  that  of  the  most  promising  NWC  control  law. 
The  latter  attempts  to  compensate  for  scale  factor  changes  by  using 
amplification  and  limiting  with  the  HCL.  Briefly,  the  strategy  here 
Is  to  use  an*>! i f ieation  to  compensate  for  low  scale  factors  and  use 
limiting  to  reduce  the  resulting  overshoot  for  high  scale  factors. 

The  responses  of  the  two  control  laws  for  both  maximum  and 
minimum  scale  factors  at  four  different  initial  attitudes  are 
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FIGURE  5.  Block  Diagram  of  HHCL  Gain  Implementation 

summarized  in  Table  3*  Also,  the  response  of  the  ideal  DCCL  to 
an  initial  roll  error  is  summarized,  Case  #  1.  This  case  was  simu- 
lated  to  provide  a  reference  response  in  the  sense  that  the  angular 
velocity  loops  were  designed  assuming  the  use  of  this  control  law, 
and  to  facilitate  selection  of  K,  Equation  (28).  Cases  #2  through 
#9  allot?  easy  comparison  of  the  HCL  with  amplification  and  limiting 
for  conditions  of  maximum  and  minimum  scale  factors.  Significant 
variations  in  the  time  required  to  come  within  3*  of  the  vertical 
are  apparent  in  eacn  pair  of  cases  (2“3»  *»“5,  6-7.  and  8-9),  '  Actual 
percentage  differences  are  even  greater,  because  tl~»s  shown  include 
approx imateiy  150  msec  of  time  spent  on  the  launch  rails,  before 
any  maneuver  could  begin.  Significant  variations  In  the  amount  of 
overshoot  are  also  evident.  The  initial  attitudes  of  Cases  #2  -  #9 
correspond  with  those  of  Cases  #10  -  #17,  which  used  the  KHCl. 
Comparison  of  each  pair  of  these  latter  cases  (10-11,  12-13,  1*»“15, 
and  16-17)  reveals  much  more  consistent  performance  than  was 
achieved  using  the  HCl  with  amplification  and  limiting.  A  compar¬ 
ison  of  the  control  laws  at  corresponding  initial  conditions  (Cases 
2  and  10,  3  and  II,  k  and  12,  5  and  13,  6  and  \k,  7  and  15,  8  and 
16,  and  S  and  17)  indicates  that  the  HHCl  Is  slower  than  the  HCl 
in  reaching  the  3*  crossover  for  most  cases.  However,  the  former 
exhibits  less  overshoot  than  the  latter  in  most  cases. 

The  important  result  of  the  above  comparison  is  not  that 
the  HCl  (with  an^)l  if  ication  and  limiting)  is  usually  faster  titan 
the  HHCL,  since  no  attempt  was  made  to  optimize  the  response  of 
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the  latter.  Rather,  the  important  result  Is  that  the  HliCL  provides 
comparable  responses  which  are  relatively  Insensitive  to  scale 
factor  changes.  It  should  be  noted  that  the  MHCL  responses  were 
obtained  with  the  estimated  output  of  an  antenna  pointing  at  the 
zenith,  in  Equation  (20),  equal  to  its  actual  value  in  the 
simulation  data.  This  approximation  should  not  significantly 
affect  the  results  as  only  small  variations  in  the  zenith  temper'* 
ature  are  anticipated.  If  the  difference  between  the  estimated 
and  actual  zenith  temperatures  were  to  become  a  problem,  it  might 
be  necessary  to  update  the  estimate  on-line. 

The  techniques  presented  in  the  section  labeled  "Additional 
Compensation"  were  not  evaluated  using  the  6-liOF  simulation. 


RECOMMENDATIONS 


The  previous  analysis  and  evaluation  not  only  establishes  the 
potential  of  the  Modified  Heuristic  Control  Law,  Equation  (20),  for 
use  with  the  present  two-axis  control  system,  but  also  suggests  cn 
increased  potential  for  dealing  with  more  optimal  maneuvers  which 
may  be  desired  if  a  three-axis  control  system  is  implemented.  As 
a  result,  it  is  recommended  that  this  control  law  be  considered  a 
leading  candidate  for  the  MPES/VSS  system.  This  recommendation  is 
made  with  the  knowledge  that  additional  effort  Is  required  to 
provide  logic  for  inverted  attitudes  and  to  optimize  performance 
for  the  range  of  flight  conditions.  It  is  also  recommended  that 
the  techniques  in  the  section  labeled  "Additional  Compensation" 
be  evaluated  for  use  with  a  three-axis  cc.itrol  system. 


FURTHER  STUDY 


Several  areas  which  appear  appropriate  for  further  study  are 
briefly  discussed. 

An  attempt  should  be  made  to  optimize  the  control  laws  for  all 
initial  flight  conditions.  This  optimization  process  should 
consider  such  things  as  probable  trajectory  of  aircraft  after 
ejection,  direction  of  relative  wind,  maximum  maneuvers  tolerated 
by  the  pilot,  and  seat  attitude  and  motion  at  the  time  of  parachute 
deployment.  Also,  the  techniques  presented  in  the  section  labeled 
"Additional  Compensation"  should  be  evaluated. 

Data  should  be  obtained  to  determine  the  MARS  output  charact- 
eristics  for  measurements  over  water  and  along  the  land-water 
interface.  These  date  should  be  used  to  re-evaluate  the  candidate, 
control  laws. 

The  number  and  orientation  of  the  MARS  antennas  should  be 
analyzed  considering  the  candidate  control  laws  to  determine  poss¬ 
ible  advantages  of  alternate  antenna  configurations. 

Current  plans  call  for  testing  and  evaluation  of  a  thrtee-axis 
control  system.  This  additional  capability  will  require  that  the 
vertical  reference  system  used  previously  be  expanded  to  an 
attitude  reference  system  which  Includes  a  heading  reference.  One 
candidate  reference  system  would  use  outputs  from  both  the  MARS  and 
strapdown  rate  gyros.  Previously,  roost  strapdown  systems  used  the 
gyro  outputs  to  update  a  3  X  3  Direction  Cosine  Matrix  (DCM)  which 
relates  the  body  axes  to  reference  axes.  Increasingly  popular 
alternatives  to  DCM  methods  Include  the  use  of  quaternions  which 
are  related  to  Euler  parameters.  These  methods  are  based  on  Euler's 
Theorem  which  states  (In  effect)  that  any  rotation  of  one  coordinate 
system  with  respect  to  another  may  be  described  by  a  rotation  of 
sonx3  angle  about  a  single  fixed  axis.  The  Euler  parameters  define 
this  rotation  angle  and  axis,  and  the  quaternion  is  a  compact  form 
for  representing  these  parameters.  Quaternions  offer  advantages 
over  the  DCH  such  as  increased  accuracy  and  decreased  computational 
requirements,  since  the  former  uses  only  four  parameters  compared 


to  nine  for  the  latter.  Also  conditions  of  nonorthogonality, 
resulting  from  computation  errors,  are  much  easier  to  correct  for 
quaternions  than  for  DCM's.  The  quaternion  Is  Ideally  suited  for 
the  attitude  control  problem  as  it  defines  a  rotation  axis  and 
angle  which  will  align  the  body  axes  with  the  reference  axes.  This 
maneuver  is  optimal  in  the  sense  that  it  minimizes  the  total  angle 
through  which  the  body  rotates.  It  is  suggested  that  quaternions 
be  considered  for  use  with  future  MPES/VSS  systems.  Quaternions 
are  discussed  further  In  Appendix  0. 
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APPENDIX  A. 

MARS  OUTPUT  DATA 

The  following  data,  provided  by  NWC,  were  obtained  by  rolling 
MARS  about  its  x  body  axis  and  recording  the  outputs  corresponding 
to  each  antenna.  Note  that  the  roll  angle  is  also  the  angle 
between  Antenna  #1  and  the  zenith. 


ROLL 

ANTENNA 

OUTPUTS 

ANGLE 

#1 

n 

#3 

#4 

•  • 

u.atOftV 

6. 439mV 

7. 013mV 

5. 992mV 

9 

1 1 . 308mV 

6. 401  mV 

6. 080m V 

6. 381mV 

It 

1 1 . 319mV 

6.909MV 

9.320MV 

7.097mV 
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11. 324mV 

6.460MV 

4*  632mV 

7.719mV 

2* 

1 1 . 337mV 

6. 477mV 

3.920MV 

6. 369mV 

29 

11 . 272mV 

6. 463mV 

3. 393mV 

8. 935mV 

30 

1 1 . 227mV 

6. 482mV 

2. 921aV 

9. 402mV 

39 

1 1.331aV 

6. 914mV 

2. 967mV 

9. 087mV 

40 

11.291 mV 

6.520MV 

2. 341 mV 

10. 220mV 

49 

1 1 . 164mV 

6. 491mV 

2.120mV 

10.470MV 

90 

10.990MV 

6.470MV 

1 . 901mV 

10. 669mV 

99 

10. 661 mV 

6.429UV 

1.838mV 

10. 796mV 

60 

10. 207mV 

6. 379mV 

1 . ?69mV 

16.901mV 

69 

9.  671  mV 

6*  373mV 

1 . ?53mV 

10. 9 99 mV 

70 

9.400MV 

6. 363mV 

1.743mV 

11 »036mV 

79 

6.81UV 

6. 393mV 

1 . 807mV 

11. 102mV 

60 

7 . 999mV 

6. 339mV 

1 . 927mV 

11.1 16mV 

89 

7.  I20mV 

6. 280m V 

2. 077m V 

II.UOmV 

98 

6.362mV 

6. 286m V 

2.2I0MV 

1 1 . 127mV 

95 

9. 931  mV 

6. 237mV 

2. 380m V 

1 1 . 094mV 

100 

4*  794mV 

6.209MV 

2. 943mV 

1 1 . 002mV 

109 

4. 094aV 

6. 179mV 

2.606mV 

1 1 . 049mV 

110 

3.  564mV 

6. 1€1mV 

2. 600m V 

1 1 . 049mV 

115 

3. 123«V 

6. 197mV 

2. 714mV 

1 1 . 824mV 

120 

2.700MV 

6. 145mV 

2  *671  mV 

10. 909mV 

129 

2.478&V 

6.094MV 

2.390mV 

10.0C0MV 

130 

2.269mV 

6. 1O0MV 

2.923mV 

10. 099m V 

139 

2. 004mV 

6.12UV 

2.952mV* 

10. 712mV 

140 

1.961 mV 

6. 119mV 

2.637mV 

10. 969mV 

149 

1 . 890mV 

6.O70MV 

2.692mV 

19. 333mV 

190 

1 . 843aV 

6.054mV 

2.929mV 

10. 1 14mV 

199 

1 . 032mV 

6. 899mV 

3.279mV 

9.757mV 

160 

1 . 013mV 

6.022mV 

1.692m V 

9. 313mV 

1 65 

1 , 779mV 

9. 969mV 

4.048MV 

0.O11MV 

170 

1 . ?39mV 

9.949mV 

4. 996m V 

8. 170mV 

179 

1 . 733mV 

5.930mV 

9. 214mV 

7.910mV 

100 

1. 721mV 

9.9I3mV 

9.801mV 

6. 032m V 
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APPENDIX  B. 

DERIVATION  OF  ANGULAR  VELOCITY  CONTROL  LOOP 
LOW  FREQUENCY  TRANSFER  FUNCTION 


A  block  diagram  of  the  control  loop  was  obtained  from  NWC  and 
reduced  to  that  shown  In  Figure  B-l  by  neglecting  the  transfer 
functions  of  the  rate  gyro,  filter  circuit,  and  part  of  the  actuator. 
The  closed-loop  transfer  function  of  this  system  Is 

SU  (S)  .  - (S  +  20) - 

q  S*  +  131  $  ♦  64l6  S  +  65320 

c 

The  roots  of  the  characteristic  equation  are  found  using  spectral 
separation.  At  high  frequencies  S  »  1  and  therefore  the  constant 
term  is  neglected  and  the  poles  are  the  roots  of 

S2  131  +  6416  -  0 
or 

S  -  -66  ±  461  ,  (B-2) 


at  low  frequencies  S  «  1  and  therefore  the  poles  are  the  roots  of 
6416  S  ♦  65320  -  0 


or 

S  =  -10  (B-3) 

The  roots  of  the  characteristic  equation  are  approximately  the  values 
given  in  Equations  (B-2)  and  (B“3)*  Since  the  roots  of  Equation 
(B-2)  are  significantly  "faster"  than  that  of  Equation  (B-3),  the 
former  were  neglected  In  obtaining  the  following  low  frequency 
transfer  function. 


S-  (5) 
q„ 


3266  (S  +  20) 
“  641 6  S  l  65320 


0.05  S  +  1 
“  0.1  S  +  1 


(B-4) 
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B-1.  Block  Diagram  of  Low  Frequency  Transfer  Functions 
in  Angular  Velocity  Control  Loop. 


(Antenna  Output)  (MV) 


APPENDIX  C. 

SCALING  OF  MARS  OUTPUT  DATA 


The  NUC  6-DOF  digital  simulation  accommodates  three  different 
sets  of  MARS  data  which  relate  the  angle  between  the  nadir  and  an 
antenna  to  the  sensed  temperature  and,  therefore,  to  an  output 
signal.  The  data  sets,  which  correspond  to  the  warmest,  nominal, 
and  coolest  expected  nadtr  temperatures,  are  obtained  by  scaling 
the  data  for  Antenna  §\  in  Appendix  A.  A  sketch  of  the  resulting 
data  is  shown  in  Figure  C-1.  Note  that  the  independent  variable 
(abscissa)  has  been  changed  to  the  angle  from  the  nadir,  In  order 
to  be  consistent  with  the  6-DOF  simulation.  The  scaling  was 
accomplished  assuming  the  data  In  Appendix  A  correspond  to  zenith 
and  nadir  temperatures  of  15*K  and  250#K>  respectively.  Also, 
the  nominal  and  minimum  nadir  temperatures  are  assumed  to  be  200*K 
and  150*K,  respectively.  The  scaled  data  set  for  the  warmest 
nadir  surface  was  obtained  by  normalizing  the  Antenna  #1  data  in 
Appendix  A  to  the  zenith  value,  11.318  MV. 


^9*®~fro£ 


FIGURE  C-1.  Sketch  of  Scaled  Output  Oata 
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Scaled  Outputs  (A) 


The  temperature  differences  between  the  zenith  and  both  the  coolest 
and  nominal  nadir  surfaces 


AT  =  150°-  15*-  135°K 

C 

AT  -  200*-  15*-  185*K 
n 


were  then  divided  by  the  assumed  difference  for  the  warmest  nadir 
surface 

AT  -  250*-  15*-  235“ K 
w 

Finally,  the  resulting  proportionality  constants  were  used  to 
scale  the  antenna  output  voltage  for  the  warmest  nadir,  in  order 
to  obtain  data  sets  for  the  nominal  and  coolest  nadirs,  as 
fol lows 

An  "  V(l8 0*7  p180*)  *  (V(180°)  -  Vw  (ajj 

and 

Ac  (a)  *  vTTBoT  (V(l80<)  *  Vw 

where 

An  (a)  -  Scaled  data  for  nominal  surface 

A  (a)  -  Scaled  data  for  coolest  surface 
c 

V  (a)  -  Antenna  output  for  warmest  nadir 
'  cond i t I on* ,  HV 

V(l80*)  -  Output  of  antenna  pointing  at  zenith, 
NV.  (Assumed  the  same  for  all  nadir 
temperatures) 

a  ■  Angle  from  antenna  axis  to  nadir 


The  resulting  data  sets  are  shown  in  Table  C-1 . 


APPENDIX  D. 

ATTITUDE  CONTROL  USING  QUATERNIONS 


One  method  of  describing  the  relationship  between  two  coordinate 
frames  is  the  3X3  Direction  Cosine  Matrix  (DCM).  Many  strapdown 
systems  use  rate  gyro  outputs  to  update  the  DCM,  thereby  maintaining 
an  attitude  reference.  An  advantage  of  this  method  is  that  the 
matrix  required  to  transform  vectors  between  the  coordinate  systems, 
that  is  the  DCM,  is  available  directly.  An  alternate  method  of 
describing  orientation  involves  the  use  of  quaternions.  This 
method  is  based  on  Euler's  Theorem  which  states  (in  effect)  that  any 
rotation  of  one  coordinate  system  with  respect  to  another  may  be 
described  by  a  rotation  a  about  a  fixed  axis  E.  In  one  form,  the 
quaternion  presents  this  rotation  as  a  four  component  vector, 
Reference  3. 

a 

COS  j 

Ex  sin-  (D-l) 

.  a 
Ey 

lz  sin  •y 


where  the  unit  vector 


The  components  of  Q  are  referred  to  as  the  Euler  symmetric  para¬ 
meters  and  were  chosen  to  simplify  computations  of  successive 
rotations,  Reference  k.  Note  that  these  components  satisfy  the 
normality  condition 

qj2  +  q22  ♦  q^2  ♦  q^2  -  1  (D~3) 
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One  method  of  combining  quaternions  to  describe  successive 
rotations  Involves  the  use  of  matrix  notation.  Consider  three 
coordinate  systems  a,  b,  and  c.  If  the  quaternion  describes 
the  rotation  of  a  Into  b,  and  the  quaternion  Qcb  describes  the 
rotation  of  b  Into  c,  then  the  rotation  of  a  Into  c  Is  described 
by  the  quaternion  (Reference  4) 


Components  of  «cb 


An  attitude  reference  can  be  maintained  by  updating  a 
quaternion  describing  the  relationship  between  Inertial  and  body 
axes  using  the  differential  equation  (Reference  5) 
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where  p,  q,  and  r  are  the  body  angular  velocities  obtained  from 
strapdown  rate  gyros. 

The  above  equations  define  and  describe  some  properties  of  a 
quaternion  represented  as  a  four  component  vector.  Some  advantages 
of  quaternions  over  the  DCH  for  use  In  attitude  reference  systems 
are  now  discussed.  First,  the  quaternion  contains  four  components 
and  Is  therefore  only  once  redundant,  whereas  the  DCH  contains 
nine  elements  and  Is  therefore  six  times  redundant.  Second,  comput¬ 
ation  errors  entering  Into  tho  DCH  require  normalization  and  ortho- 
gonalization  of  three  vectors.  The  latter  process  can  be  difficult. 
The  corresponding  errors  in  a  quaternion  merely  require  normal  Ization 


to  satisfy  the  condition  of  Equation  (0-3) •  Third,  successive 
rotations  using  DCH's  require  multiplication  of  two  3X3  matrices, 
which  leads  to  27  multiplications  and  18  additions,  compared  with 
16  multiplications  and  12  additions  for  the  same  operation  using 
quaternions,  Equation  (0-4).  Fourth,  the  strapdown  differential 
equation  tor  a  DCM  requires  18  multiplications  and  9  additions, 
compared  with  12  and  8  respectively  for  quaternions,  Equation  (D-5)» 
For  systems  which  require  the  transformation  of  vectors  between 
inertial  and  body  axes,  the  above  computational  advantages  of  quat¬ 
ernions  may  be  lost  when  the  OCH  is  calculated  from  the  quaternion. 
However,  most  attitude  control  systems  do  not  require  direct  vector 
transformations  and,  therefore,  do  not  require  explicit  knowledge 
of  the  DCM.  In  these  cases  the  quaternion  retains  its  computational 
advantages. 

Another  advantage  of  quaternions  is  that  they  explicitly  define 
a  single  rotation  axis  and  angle  which  will  result  in  alignment  of 
the  inertial  and  body  coordinate  systems.  If  the  angular  velocity 
vector  Is  along  this  axis,  the  resulting  maneuver  will  be  optimum 
in  the  minimum  angle  sense.  Therefore,  a  logical  approach  Is  to 
define  angular  velocity  commands  in  terms  of  the  quaternion  para¬ 
meters.  Possibilities  include 
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where  the  subscript  c  indicates  commanded  values  and  K  is  a  constant 
gain. 

Initialization  of  the  quaternion  or  OCK  is  required  if  either 
is  used  in  a  strapdown  attitude  reference/control  system.  For 
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3-axis  control  of  the  MPES/VSS  escape  system  it  Is  desired  to  use 
the  MARS  sensor  to  perform  this  Initialization.  This  approach  is 
complicated  by  the  limitations  of  the  MARS  outputs.  Previous 
analysis  has  shown  that  MARS  can  provide  a  reference  to  the  vertical, 
however,  it  has  not  indicated  how  to  extract  (changes  in)  yaw 
attitude.  Therefore,  yaw  attitude  must  be  referenced  to  its  initial 
value,  and  changes  from  this  value  should  either  be  calculated  using 
the  rate  gyro  outputs  or  minimized  by  proper  selection  of  the 
commanded  angular  velocity.  One  method  of  accomplishing  the  latter 
is  to  estimate  the  pitch  and  roll  angles  using  the  MARS  outputs, 
form  quaternions  for  first  returning  the  roll  attitude  to  zero 
and  then  returning  the  pitch  attitude  to  2ero,  and  use  Equation 
(0-4)  to  calculate  the  quaternion  for  these  successive  rotations. 

This  resulting  quaternion  defines  a  single  axis  rotation  which 
will  align  the  seat  with  the  vertical  without  changing  its  heading. 
This  initialized  quaternion  can  be  propagated  by  numerical  inte¬ 
gration  of  Equation  (0-5)  and  a  control  law  similar  to  those  in 
Equations  (0-6)  and  (0-7)  can  then  be  used.  A  second  method  is 
to  estimate  the  minimum  rotation  axis  and  angle  to  the  vertical 
using  the  MARS  outputs,  form  the  corresponding  quaternion  and 
determine  the  associated  heading  change,  form  a  quaternion  to 
return  the  heading  to  its  original  value,  and  calculate  the  quat¬ 
ernion  for  these  successive  rotations.  This  quaternion  should 
correspond  to  that  obtained  using  the  previous  method  and,  therefore 
can  be  used  in  the  same  manner. 

in  the  approaches  discussed  above  the  KARS  outputs  were  only 
used  to  initialize  a  quaternion.  An  alternate  approach  is  to  also 
use  the  MARS  outputs  to  update  the  attitude  estimate  at  given  time 
intervals  along  the  trajectory.  A  standard  approach  to  this  problem 
of  combining  attitude  and  angular  velocity  data  is  to  use  a  steady- 
state  Kalman  filter. 

For  some  flight  conditions,  limited  control  authority  may  make 
it  undesirable  to  command  an  angular  velocity  along  the  "Quaternion 
axis."  in  these  cases  the  original  quaternion  can  be  propagated 
and  a  modified  quaternion  used  for  control.  At  a  later  point  along 


tl s  trajectory,  control  can  be  switched  to  the  original  quaternion 
which  will  provide  an  optimum  maneuver  to  the  vertical  and  original 
heading.  The  modified  quaternion,  mentioned  above,  may  result  from 
simple  changes  to  the  original  or  may  result  from  a  successive 
rotation  calculation,  depending  on  the  desired  control  strategy. 

The  above  methods  of  using  the,  MARS  outputs  to  initialize 
quaternions  should  be  evaluated  with  respect  to  their  relative 
accuracies  and  ease  of  implementation  in  a  microprocessor.  Also, 
advantages  of  using  the  MARS  outputs  to  update  the  seat  attitude 
should  be  evaluated.  Finally,  modified  quaternion  control  laws 
should  be  formulated  and  their  effectiveness  In  dealing  with  limited 
control  authority  evaluated.  These  evaluations  should  include  the 
use  of  suitable  simplifying  approximations. 
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